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Abstract. We present sub-milliarcsecond observations of SiO masers in the late-type stars IRC +10011 and \ Cyg- 
We have used the NRAO Very Long Baseline Array (VLBA) to map the 43 GHz («=1,2 J=l-0) and the 86 
GHz (u=l,2 J=2-l) SiO masers. All the transitions have been imaged except the v=2 3=1— \ in IRC +10011. 
We report the first VLBI map of the v=l J=2-l 28 SiO maser in IRC +10011 as well as the first VLBA images 
of SiO masers in an S-type Mira variable, \ Cyg- In this paper we have focused on the study of the relative 
spatial distribution of the different observed lines. We have found that in some cases the observational results 
are not reproduced by the current theoretical pumping models, either radiative or collisional. In particular, for 
IRC +10011, the v—1 J=l-0 and J— 2-1 28 SiO lines have different spatial distributions and emitting region sizes, 
the J=2-l emission being located in an outer region of the envelope. For y Cyg, the distributions also differ, 
but the sizes of the masing regions are comparable. We suggest that the line overlaps between ro-vibrational 
transitions of two abundant molecular species, H2O and 28 SiO, is a possible explanation for the discrepancies 
found between the observations and the theoretical predictions. We have introduced this overlapping process in 
the calculations of the excitation of the SiO molecule. We conclude that the line overlaps can strongly affect the 
excitation of SiO and may reproduce the unexpected observational results for the two sources studied. 

Key words, masers -technique: interferometric - stars: circumstellar matter- stars: AGB and post-AGB 



^ 1. Introduction 



Circumstellar SiO maser emission is detected in Long 
Period Variables (AGB-stars): Mira-type, semi-regular 
and irregular variables, and OH/IR stars; of M and S spec- 
tral types. Due to the strong energy requirements for the 
pumping of these masers, they are produced in the in- 
nermost shells of the circumstellar envelopes around these 
stars, where the gas density and stellar radiation field are 
important. On the other hand, this amplified emission is 
very bright and compact, therefore being an ideal target 
for Very Long Baseline Interferometric (VLBI) studies. In 
fact, VLBI observations are the only way to accurately 
map the structure of the SiO maser emitting regions, pro- 
viding angular resolutions up to several 0.1 mas (Colomer 
et al. nMty . 

First estimations of the location of the SiO maser 
emission were derived from single-dish observations and 
from theoretical modelling (see e.g. Alcolea et al. 
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Bujarrabal et al. I1981f) . These works already suggested 
that the SiO masers originate at distances of a few stellar 
radii from the central star, between the extended atmo- 
sphere and the dust-formation point. This has been con- 
firmed in subsequent VLBI measurements (Danchi et al. 
ITggi Greenhill et al. \HM\ Wittkowski & Boboltz E^tMjl 
and in addition, that the distribution of the emission con- 
sists of ring-like structures composed of multiple maser 
spots with typical sizes of a few mas (Diamond et al. Q.994 
Greenhill et al. ITM51 Desmurs et al. I2UUU|) . VLBI moni- 
toring of the SiO maser emission probes the kinematics 
of the gas in these inner shells, revealing inward/outward 
motions (Boboltz et al. lTW7l Diamond & Kemball I20"0"3|l 
or showing indications of rotation of the maser shell in 
some peculiar objects (Boboltz & Marvel 2000; Hollis et 
al-EOOIl Sanchez Contreras et al.EEEl- 

These observations also provide very valuable con- 
straints for the theoretical models accounting for the ori- 
gin and pumping of these masers, which are needed in or- 
der to understand the mechanism itself and to derive the 
physical conditions of the maser emitting regions from the 
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Table 1. Properties of the observed sources 



source 


R.A. 


Dec. 


variability 




period b 


phase c 


dist. d 


^LSR E 




(J2000) 


(J2000) 


type 


(^yr- 1 ) 


(days) 




(kpc) 


(kms^ 1 ) 


IRC +10011 


01 h 06 m 25 ! :99 


+12°35'53'.'3 


OH/IR 


8.1 10" 6 


630 


0.10 


0.47 


9.2 


xCyg 


19 h 50 m 3S92 


+32°54'50'.'6 


Mira 


4.1 10~ 7 


408 


0.25 


0.11 


10.2 



a Mass loss values derived following Loup et al. (11993t but for the distances adopted here. 
b Periods form Pardo et al. (I2UU4II . 

c Stellar phases at the time of the observation: for IRC +10011, extrapolated value from Pardo et al. |2004| with respect to 
the NIR maximum (the typical phase lag between IR/SiO and optical maxima is ~0.1 periods). For xCyg, the optical phase 
derived from the light curve provided by the AAVSO (see http://www.aavso.org). 

d Distances from Herman et al. (THjSBI for IRC +10011 and from the HIPPARCOS catalogue for xCyg. 

c Systemic velocity from high resolution CO observations by Cernicharo et al. 1199711 and Bujarrabal et al. 119898 respectively. 



observations. These theoretical models are divided in ra- 
diative and collisional, the main differences between both 
types being the main process that provides the energy 
to maintain the inversion of the population and therefore 
produce the maser transitions. The responsible mechanism 
in the case of the radiative models is the IR radiation from 
the central star (e.g. Bujarrabal 19!) la 10!) lb), whereas in 
the collisional models the energy is obtained from the ki- 
netic energy of the gas (e.g. Lockett & Elitzur 1992). 

Most of the VLBI observational results on SiO masers 
have been obtained mainly in M-type Mira stars (sim- 
ply because they typically present stronger maser emis- 
sion than the other types of variables) and from observa- 
tions of the 7 mm rotational transitions of 28 SiO: v=l and 

2 J— 1-0 lines. Successful VLBI observations of the 3 mm 
J=2-l rotational transition were made by Doeleman et al. 
ljT55£)l and Phillips & Boboltz pitH)) towards the late type 
stars VX Sgr and o Cet, which helped in the study of the 
spatial structure of these 3 mm SiO masers. Nowadays, us- 
ing the NRAO 1 Very Long Baseline Array (VLBA), which 
was upgraded to operate at 3 mm wavelength in 2001, it 
is possible to compare the observed results at 3 mm with 
those at 7 mm (Phillips et al. I2003|) . As we will see, the 
comparison of these quasi simultaneous observations with 
the theoretical predictions for the J=2-l lines will fur- 
ther constrain the models proposed to explain the physics 
involved in the SiO maser phenomenon. 

Making use of its new observational capabilities at 

3 mm, we have started multi-epoch/multi-transitional 
high resolution observations of SiO masers using the 
VLBA in a sample of long period variable stars. We 
have studied the SiO maser emission in the sources 
IRC +10011, xCyg and TXCam. In particular we have 
observed the v=l, v=2 J— 1-0 and J=2— 1 transitions of 
28 SiO and the v=0 J=l-0 and J=2-l lines of the 29 SiO 
isotopic substitution (see Sect. 2). 

We present in this paper the results obtained in the 
first epoch for IRC +10011 and xCyg (those for TXCam 

1 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative 
agreement by Associated Universities, Inc. 



will be presented in a separate paper). We report the first 
VLBI map of the v=l J =2-1 28 SiO maser in IRC +10011 
as well as the first VLBA images of the v—1 and v=2 
J =1-0 and J=2— 1 transitions in an S-type Mira variable, 
X Cyg. We have detected all the observed lines except the 
28 SiO v=2 J =1-1 in IRC +10011 and, unfortunately, the 
29 SiO lines in the two variables. A detailed description of 
the different maps obtained is given in Section 3. 

In Section 4, we study and compare in detail the dif- 
ferences between the maser emissions at 7 mm (v=l and 
v=2 J=l-0) and between the 7 and 3 mm (v=l and 2, 
J =2—1 and J=l-0). We focus on their relative spatial 
distribution, in order to compare the observational results 
with the predictions of the currently proposed pumping 
models for circumstellar SiO masers. We have found that 
our results for the 3 and 7 mm comparison are not repro- 
duced by any of these models, since all predict a similar 
emitting region at both wavelengths (Bujarrabal 1994a 
Humphreys et al. .20QDI and we however obtain a totally 
different distribution in both IRC +10011 and xCyg. We 
propose the line overlaps between ro-vibrational transi- 
tions of two abundant molecular species, H2O and 28 SiO, 
as a possible explanation of those unexpected results. For 
this reason, we have also introduced this overlapping pro- 
cess in the calculations and discuss the physical effects. 
Finally, the main conclusions and results derived from the 
high resolution observations and theoretical calculations 
performed are summarized in Section 5. 

2. Observations and data analysis 

We have performed sub-milliarcsecond resolution observa- 
tions of SiO masers transitions with the VLBA. We have 
observed the long period variable stars IRC +10011 on 
2002 January 5 and xCyg on 2001 May 9 (Table [TJ. We 
have measured the 28 SiO rotational lines J=l-0 and 3=2- 
1 in the v=l and v=2 vibrationally states and the J=l— 
and J =2-1 transitions of the 29 SiO isotopic substitution 
in the ground vibrational state v=0. 

All available antennas were used at the time of the 
observations, the entire array (10 antennas) for the 7 mm 
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Table 2. Observational setups 



setup 


number 
of IFs 


BW 
(MHz) 


spectral 
channels 


spectral 
resolution (km/s) 


SiO transitions observed 
(Left circular polarization only, except when noted) 


7 mm 


4 


8 


256 


0.22 


28 SiO v=l J=l-0 t ; 28 SiO v=2 J=l-0; 29 SiO v=0 J=l-0 


3 mm-A 


2 


16 


512 


0.11 


28 SiO v=2 J=2-l; 29 SiO v=0 J=2-l 


3 mm-B 


2 


16 


512 


0.11 


28 SiO v=l J=2-l f 



Left and right circular polarizations (LCP & RCP) . 



Table 3. Summary of the observational results 





IRC +10011 


xCyg 


band species maser rest free. 

transition (MHz) 


restoring # of masing region 
beam (mas) spots size (mas) 


restoring # of masing region 
beam (mas) spots size (mas) 


7mm 28 SiO v=l J=l-0 43122.080 
v=2 J=l-0 42820.587 
29 SiO w=0 J=l-0 42879.916 


1.1x0.3 9 19-24 
0.4x0.2 11 16-21 
non-det. 


0.7x0.2 6 50-62 
1.0x0.6 1 
non-det. 


3mm 28 SiO c=l J=2-l 86243.442 
v=2 J=2-l 85640.456 
29 SiO w=0 J=2-l 85759.132 


0.7x0.3 8 28-33 

non-det. 

non-det. 


0.7x0.5 6 50-58 

0.7x0.3 4 

non-det. 



(43 GHz) and 7 antennas for the 3 mm (86 GHz) obser- 
vations. The distances between the VLBA antennas range 
between 236km (LA-PT), and 8611km (SC-MK) for the 
7 mm and 6156 km (NL-MK) for the 3 mm observations, 
thus achieving maximum angular resolutions of about 0.17 
and 0.12 mas respectively. 

The data were correlated at the VLBA correlator in 
Socorro (New Mexico, USA). The observations for each 
of the sources consisted of 3 different lines setups, the 
7 mm, the 3 mm-A and 3 mm-B, which comprised a total 
observing time of about 8 h per source. A summary of the 
observing and correlator setups is presented in Table |3 
The u-v coverages are shown in Fig.^ 

The observations carried out are summarized in 
Table All the 28 SiO transitions observed were detected 
except the v=2 J=2-l maser in IRC +10011 and, unfor- 
tunately, none of the 29 SiO lines because their flux was 
below the detection limit (they were not detected in the 
auto-correlations either) . 

Using the mentioned setups, the four 7 mm masers 
were observed simultaneously and the 7 and 3 mm sep- 
arated only by a few hours. This is important since one of 
the characteristics of SiO masers is their time variability. 
From monitoring programs of Mira and other types of reg- 
ular variables, two types of variability have been reported: 
long-term variability on time scales of several months (see 
Alcolea et al. 119991 and references therein), and changes 
on shorter time scales, 10-20 days (Pijpers et al. 11994(1 . 
So far no significant variability has been observed on time 
scales of hours, and therefore our maps of the 7 and 3 mm 
masers can be directly compared. 



The data reduction was performed using the 
Astronomical Image Processing System (AIPS) package. 
The calibrator used to calculate the bandpass and single- 
band delay corrections, which was later applied to the 
spectral line sources, was 3C 454.3 for both sources. 

The amplitude calibration was done by comparing a 
calibrated total power spectrum of the line source, used 
as a template, with the autocorrelations measured at each 
antenna, except in the v—1 J=2-l maser in IRC +10011 
and in the v—1 J=2—l transition in x Cyg. In these cases, 
we used the system temperatures and antenna gains tables 
provided by the NRAO within the data files. 

The determination of the phase calibration was done 
following the standard method for spectral line VLBI 
data. Residual phase errors were removed later in a self- 
calibration process. Finally, the integrated intensity maps 
were obtained independently for each of the SiO transi- 
tions detected, using the image deconvolution CLEAN al- 
gorithm. 

Note that when the phase calibration is done by means 
of phase-referencing techniques, the absolute positions in 
the maps are lost and only relative positions to the ref- 
erence channel are obtained. This adds a difficulty when 
comparing two maps since their relative position is un- 
known. 

The results obtained are presented in Figs. 2-7, where 
for each of the mapped SiO transitions, we show: 1) The 
VLBA integrated intensity map (in Jy per clean beam). 
2) The individual spectra of all the maser features with 
a signal-to-noise ratio larger than 5 and appearing on at 
least three consecutive spectral channels. The spectrum 
for each spot was calculated transforming the map units 
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(Jy/beam) to flux units (Janskys) for the selected maser 
region. 3) The total power calibrated spectrum of the 
maser transition. 4) The total flux recovered in each fre- 
quency channel after the calibration and imaging. 5) The 
ratio of the total power to the recovered flux. 

3. Results 

In this section, we briefly describe the most relevant as- 
pects of the observational results obtained for the 28 SiO 
masers observed in the long-period variables IRC +10011 
and x Cyg. 

We note that there is an additional difficulty when esti- 
mating the angular extents of the different masing regions 
since, in some cases, the distributions are not complete 
rings. For this reason, we give an approximate range of 
the angular sizes derived from the brightest components 
that appear in the maps. 

3.1. IRC +10011 

IRC +10011 (WX Psc) is an OH/IR star of spectral type 
M9, characterized by a very thick circumstellar envelope 
and higher mass loss rate than Mira-type variable stars 
(see TableHJ. The SiO emission of IRC +10011 varies in 
time with a period of about 630 days (Herman & Habing 
1985; Pardo et al. l2004J) : our observations were performed 
close to a maximum. All the 28 SiO masers observed were 
detected except the v—2 J— -2-1, probably due to the weak 
emission of this line in oxygen-rich variable stars (see Sect. 
4.5). 

3.1.1. The v=l J— 1-0 transition 

The v=l J=l SiO results are displayed in Fig.[2 The 
intensity maps of the LCP and RCP observations were 
independently produced, although the same template au- 
tocorrelation spectra and phase reference channel (corre- 
sponding to a velocity of 9.16kms _1 ) were used in both 
cases to calibrate the amplitude and the residual fringe- 
rates. No difference was found between the two circular 
polarizations and so the two maps were finally averaged. 

The interferometric VLBA map shows nine maser 
spots distributed in three different regions with velocities 
ranging from 6 to 12kms _1 and peak intensities up to 9 
Jy (spot 7). The angular diameter of the emitting region is 
~ 19-24 mas, derived from the brightest features. Almost 
40% of the flux has been recovered in the interferometric 
maps at velocities near the systemic, 9.2 km s . 

3.1.2. The v=2 J=l-0 transition 

The v—2 J=l-0 results for IRC +10011 are presented in 
Fig-El The intensity map is composed of eleven maser 
clumps distributed into three different regions. The ve- 
locity range varies from 6 to 12kms -1 and the intensities 
from ~ 1 Jy (spot 4) to ~ 12 Jy (spot 9, where the in- 
tensity has been divided by two in the plot to ease the 



comparison with the rest of the spectra). For this transi- 
tion, the diameter of the emission is about 16-21 mas. At 
velocities around 7kms~ 1 , about 40% of the maser flux 
has been recovered. 

3.1.3. The v=l J=2-l transition 

The v=l J=2-l SiO line is presented in Fig. 01 We found 
eight different maser spots, with velocities from 6 to 
10.5 kms -1 and intensities up to 11 Jy (spot number 9). 
In this transition, about 13% of the emission has been re- 
covered at velocities ^7 kms -1 . The maser components 
are distributed in two regions with an angular separation 
of about 28-33 mas. 

3.2. X Cyg 

X Cyg is a Mira variable star with a mass loss rate similar 
to that typically found in other regular variables, and a 
period of about 408 days. This object was also observed 
near an SiO maximum (tp — 0.1-0.15, Alcolea et al. ^999); 
however, we note that in this case the variability of the 
SiO maser emission is not as regular as in other Mira 
variables (Pardo et al. |2l,)l) 1 i . x^yg is an S7-9 spec- 
tral type star, thus characterized by a peculiar chemical 
composition. S-type variables have similar Carbon and 
Oxygen abundances, resulting in a lower production of 
O-bearing molecules other than CO. In particular, for 
xCyg, the C/O ratio is found to be ~0.95 (Keenan & 
Boeshaar 198QJ, a value relatively low for an S-type star. 
Nevertheless, in this source the H2O abundance might be 
significantly lower compared to O-rich envelopes, since no 
H2O maser emission has ever been detected. This aspect 
is relevant when explaining the pumping of SiO masers, 
since, as we will see in Section 4.5, the overlap between in- 
frared lines of SiO and H2O strongly affects the excitation 
of the former species. 

3.2.1. The v=l J=l-0 transition 

The v=l J=l— SiO maser for x Cyg is presented in Fig.0 
Both, LCP and RCP observational data were averaged, af- 
ter the calibration to produce the final map, since no sig- 
nificant differences were found between the two polariza- 
tions. Six maser spots appear in the map, with velocities 
from 9.8 kms -1 to 14kms _1 , distributed in two regions 
separated about 50-62 mas. The brightest feature (spot 1) 
has an intensity of about 8 Jy. All the flux has been recov- 
ered at velocities near 15kms _1 , and an average of 50% 
for the other velocity channels. 

3.2.2. The v=2 J=l-0 transition 

The v=2 J=l-0 transition is shown in Fig.HJ Only one 
spot was found in the map, which is consistent with the 
simple line profile. Its velocity spread is 12-14 kms" 1 , and 
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Fig. 1. u-v coverage of the VLB A observations for IRC - 
frequency setups. 



10011 (upper panels) and x Cyg (lower panels) for the different 





7 1SB (km s l ) 
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Fig. 2. The v—1 J=l-0 28 SiO transition for IRC +10011. The figure consists of the integrated intensity map in Jy 
beam -1 - kms -1 units (center panel), the spectrum of the maser spots that appear on the map (numbered boxes), 
the total (AC) and recovered flux (XC) spectra (left panel above the map), and their ratio (right panel above the 
map). In all the spectra we indicate the systemic velocity with a dashed line. The peak intensity of the map is 7.21 Jy 
beam" 1 -kms" 1 and the contour levels are 0.26, 0.44, 0.62, 0.80, 0.98, 1.16, 1.34, 1.52, 3.04 and 6.08 Jy beam -1 • kms" 1 . 



its intensity is up to 
been recovered. 



Jy. Almost 50% of the emission has 4. Comparison of the different SiO transitions 



3.2.3. The v=l J=2-l transition 

The v=l J=2-l SiO maser is presented in Fig.[7| The 
VLBI map has six maser spots delineating a ring with 
a diameter of 50-58 mas, and with velocities from 10 to 
15 kms -1 . The brightest spot has a flux density of about 
25 Jy. Only 20% of the emission was recovered in most 
velocity channels. 



3.2.4. The v=2 J=2-l transition 

The first VLBA map of the v=2 J=2-l 28 SiO transition 
is presented in Fig-EI The image is composed of four spots 
with velocities ranging from 12 to 17 kms" 1 and flux in- 
tensities up to 2.5 Jy. For this maser, more than 95% of 
the emission is lost in the interfcrometric map. 



As noted previously, most of the properties of circumstel- 
lar SiO masers derived from VLBI studies have been ob- 
tained on M-type Mira variables and, in particular, from 
the observations of the 7 mm masers, mainly the v=l 
J=l— 0. However, one of the peculiarities of these stellar 
SiO masers is that they occur in a variety of lines, with 
different frequencies and excitation energies. Therefore, 
the comparison of the simultaneous observations harbors 
a crucial piece of information to understand the SiO maser 
phenomenon as well as to better constrain the theoretical 
models that try to explain the pumping of these masers. 
To date, most of the VLBI multi-transitional studies of 
SiO masers have focused on the comparison between v=l 
and v=2 J=l-0 lines (at 7 mm), while the comparison of 
SiO masers at different frequency bands remains almost 
unexplored. 

In this section we compare the distribution of the ob- 
served 28 SiO maser lines in IRC +10011 and xCyg. We 
discuss the differences found between our observational 
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results and what is predicted by the theoretical models, 
which are found to fail mainly in reproducing the differ- 
ences observed between the 7 and 3 mm masers. Finally, 
we explore the influence of the line overlaps between ro- 
vibrational transitions of H2O and 28 SiO in the excitation, 
and localization in the envelope, of the SiO masers. 



4.1. The v=l and v—2 J=l-0 masers 
4.1.1. Previous results 

The first comparison of VLBI observations of the v=l and 
v=2 28 SiO masers, using the KNIFE interferometer, was 
performed by Miyoshi et al. I|1994|l . concluding that these 
two lines were spatially coincident. More accurate VLB A 
observations have demonstrated that there is a separa- 
tion between the v=l and v=2 ring structures of several 
milli-arcsec, with the v=2 maser line in most cases located 
inside the v=l ring (Desmurs et al. 19991 12007)1 see also Yi 
et al. I2()()()j> . Recently, Cotton et al. I)2()()4|l have further- 
more demonstrated that this relative location of the two 
7 mm maser lines is independent of the optical phase of 
the star. We have found in our simultaneous observations 
a similar result for IRC +10011, but not for %Cyg as we 
will see in detail in the forthcoming subsections. 



4.1.2. IRC+10011 

The v=l and v=2 J=l-0 SiO maps for IRC+10011 are 
shown in Figs .[3 and Clearly, the spatial distribution 
of the two 7 mm SiO lines is quite similar, in both cases 
the maser spots being grouped in three different regions. 
However, the emission of the v—2 line appears to be dis- 
tributed in a region ~ 10% smaller than the v=l one, 
16-21 mas vs. 19-24 mas, thus arising from a slightly in- 
ner region. 

As we have mentioned in Section 2, since the absolute 
position of the maps is lost in the calibration process, their 
relative alignment is arbitrary. In this case, although as- 
suming that the coincidence of the centroids of the two 
maps seems a reasonable choice, there is an uncertainty 
in the centroid position of 2-3 mas. In Fig.OHwe present a 
possible alignment based on the fact that the spot 9 of the 
v=l map shows a velocity profile in between spots 9-11 
of the v=2 line. If this is the right alignment, one should 
conclude that there are no spatially coincident spots be- 
tween the two 7 mm masers, in agreement with Desmurs 
et al. I|2000j> . However, we must admit that other rela- 
tive positionings are possible. Nevertheless, taking into 
account that it is very unlikely to have spatially coinci- 
dent spots with different velocity profiles, we can only 
match 2-3 spots at most (out of ~ 10 detected in both 
maps). Therefore, at least 70% of the spots found in one 
transition do not appear in the other. 
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Fig. 4. Same as Fig.[2]for the v=l J=2-l SiO transition in IRC +10011. The peak intensity of the map is 5.71 Jy 
beam^-kms" 1 and the contour levels are 0.12 (only in spots 1 to 3), 0.20, 0.28, 0.36, 0.44, 0.88, 1.76 and 3.52 Jy 
beam • kms -1 . 



4.1.3. X Cyg 

As noted in Section 3.2.2, only one feature was found in 
the v=2 J =1-0 line, and so the comparison with the v=l 
is not relevant for \ Cyg, since we can not make any rea- 
sonable assumption of the relative position between the 
two maps. From the velocity profiles, the spot detected 
in the v=2 line could be spatially coincident with spots 
either 2 or 4 in the v=l, or with none of them (see Figs.0 
and^J) . Nevertheless, since we have detected 6 spots in the 
v=l J=l-0 maser, we should conclude that at least ~ 80% 
of the regions emitting in the v=l have no counterpart in 
the v=2 transition. A similar conclusion can be reached 
just by comparing the total flux profiles. 

4.2. The v=l and v=2 masers: J =1-0 vs. J =2-1 
4.2.1. Previous results 

The only published comparison of quasi simultaneous 
VLBI observations of 7 and 3 mm SiO maser lines has been 
performed by Phillips et al. (I2003|) in R Cas, which is an 
O-rich Mira-type variable. This source has a period of 430 
days, and a mass loss rate of~1.110~ 6 Mq yr _1 . These 
authors observed the v=l J=l-0 and J=2— 1 masers of 
28 SiO, the former line being stronger and more complex 
at that epoch. More than ten J=l-0 spots delineate an 
incomplete circle. The 3=2-1 line map consisted only of 



three spots, being two components possibly spatially co- 
incident in both transitions. 

4.2.2. IRC +10011 

The maps obtained for the v=l J=l-0 and J=2-l in this 
source are displayed in Figs .[21 and 0] As it can be easily 
seen, the spatial distributions of both maser emissions are 
totally different. The separation between the two groups 
of spots in the 3 mm line is much larger than the diameter 
of the maser emission ring in the J=l-0 line, 28-33 mas 
vs. 19-24 mas respectively. Assuming that the centroid of 
the emission is the same in both cases (see Fig. llQI) , we can 
conclude that the J=2— 1 originates from a shell located 
50% further away than the J=l— 0, having no spots in 
common. This is not surprising since the line profiles arc 
not very similar (contrary to what happens between the 
two 7mm lines in this source). In fact, choosing a different 
alignment for the maps, it is impossible to obtain more 
than one coincident spot. 

4.2.3. X Cyg 

When comparing the maps for the two v=l, J=l-0 and 
J =2—1, in xOyg (see Figs .[5] and , we found a result 
somewhat different than that for IRC +10011. Although 
the spots distribution is very different in the two transi- 
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Fig. 5. Same as Fig.0 for the v=l J=l-0 SiO transition in xCyg. The peak intensity of the map is 3.03 Jy 
beam -1 - km s -1 and the contour levels are 0.29, 0.39, 0.49, 0.68, 0.88, 1.08, 1.27 and 2.54 Jy beam -1 • km s -1 . 
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Fig. 6. Same as Fig.[2] for the v—2 J— 1-0 SiO transition in xCyg. The peak intensity of the map is 3.19 Jy 
beam -1 -kms -1 and the contour levels are 0.59, 0.98, 1.37, 1.76 and 2.15 Jy beam -1 - kms -1 . 
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Fig. 7. Same 
beam • kms 



as Fig.|21 for the v=l J=2— 1 SiO transition in x^yg. The peak intensity of the map is 18.36 Jy 
-1 and the contour levels are 3.80, 5.43, 7.06, 8.69, 11.95 and 13.58 Jy beam -1 - kms -1 . 
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Fig. 8. Same as Fig.0 for the v=2 J=2-l SiO transition in xCyg. The peak intensity of the map is 1.35 Jy 
beam -1 - kms -1 and the contour levels are 0.29, 0.33, 0.37, 0.43, 0.52 and 1.05 Jy beam -1 - kms -1 . 
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Fig. 9. Comparison of the v=l (greys) and v~2 (continuous contours) J=l transitions for IRC +10011. The zoom 
in plots are also displayed for the three masing regions. In the (0,0) region it has been represented the velocity averaged 
image (first panel) and the individual spectral emission channels covering the velocity range 10.9-7.2kms _1 . 
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Fig. 10. Comparison of the v—1 J— 1-0 (greys) and J=2-l (continuous contours) transitions for IRC +10011. 



tions, the size of the emitting regions is comparable. This 
can be more easily appreciated in Fig. llll in which we have 
again assumed that the centroid of the two emissions is the 
same. 

The VLB A map of the v=2 J— 1-0 maser is shown 
in Fig.|5j We also present the first VLBI map of the 
v=2 J=2-l 28 SiO maser in an AGB star in Fig.EI The 
v~2 J=2—l maser emission is confined into an elongated 
region of 7 mas, displaying a clear velocity gradient of 
3.5 km s -1 (~0.5 km s^" 1 per mas). These elongated struc- 
tures are typical in other SiO maser sources (see the case 
of TXCam, Diamond & Kemball lM^ . Due to this small 
spatial extent, it is very difficult to conclude a reason- 
able alignment for this map with respect to the v=l line 
maps. This is also the case of the other v=2 transition, the 
,7=1-0 line, which consists of a single spot. Nevertheless, 
the comparison of the two emitting regions is not obvious 
since for the ,7=1-0 there is only one maser component. 
However, its velocity distribution is quite similar to that 
of the spot number 4 of the J— 2-1, and so both emissions 
could arise from the same region of the circumstellar enve- 
lope. However, where this region is located in the envelope 
is very difficult to say, since the observed velocity range in 
the v=2 lines, ~ 12-16 km s _1 , is shared by several spots 
of the v=l masers. 

As we see, the available data to compare the 7 and 
3 mm masers are very scarce: this comparison has only 
been done in three sources, the two presented here, 
IRC +10011 and xCyg, and RCas by Phillips et al. 



(2003). However, the results obtained can not be more dis- 
parate. For the two "normal" stars, i.e. O-rich, we have one 
case, IRC +10011, in which the J =2-1 and J=l-0 v=l 
masers arise from two totally different regions, with proba- 
bly different physical conditions, whereas in RCas Phillips 
et al. ( 2003 ) claim that the masers originate from the same 
layer. For \ Cyg, which on the contrary is an S-type star, 
the situation is more similar to RCas. Nevertheless, we 
should note here that, as in the v=l and v=2 masers, 
from the few existing observations, it seems that it is diffi- 
cult to find maser spots emitting in both v=l J=l—0 and 
J=2-l lines. 



4.3. Comparison with theoretical models 

In an attempt to explain the existence and mainte- 
nance of SiO masers in the circumstellar envelopes of 
long period variable stars, several theoretical models have 
been proposed. Currently, two types of pumping mecha- 
nisms are considered, which are termed radiative (Kwan 
& Scoville H371 Deg uchi & Ig uchi 1X5751 Buj arrabal & 
Nguyen-Q-RieuHHED Buiarrabal fl994al IT991E1 Rausch et 
al. IT555j) or collisional (Elitzur 115801 Lockett & Elitzur 
IT5521 Doel et al. 115551 Humphreys et al. IT555I l2T?02|l . de- 
pending on which is the primary mechanism responsible 
for this amplified emission. For radiative models, the en- 
ergy of the maser pumping comes from the stellar radia- 
tion at 8 /mi, which corresponds to Av=l ro- vibrational 
SiO transitions. In collisional models, the kinetic energy 
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Fig. 11. Comparison of the v=l J=l-0 (greys) and J=2-l (continuous contours) transitions for xCyg. 



of the gas pumps the masers by means of collisions with 
the H2 molecules. 

In all cases, the inversion of the levels in the maser 
transition is due to the same effect, the self-trapping of 
the photons in the Av=l ro- vibrational transitions when 
these lines are opaque (see Kwan & Scoville Il974|l . The 
main de-excitation path for molecules of SiO in the vibra- 
tional excited estates (v=l and above) is the spontaneous 
decay via Av=l transitions (v—l-Q, v=2-l, etc.). The cor- 
responding de-excitation rate is proportional to 

A„/f3^+A + /f3+ 

where A- and A + are the Einstein's A-coefhcients for the 
spontaneous emission in the R- and P-branches, and (3- 
and (3 + are the corresponding escape probabilities. When 
the opacity is large, the escape probability in the envelope 
is given by the inverse of the opacity of the corresponding 
transition, which itself is proportional to the J value of the 
upper state. In these conditions, the higher the J value the 
more difficult for an SiO molecule to de-excite via Au=l 
ro- vibrational transitions. This results in a population in- 
version along the rotational ladder of the corresponding 
vibrationally excited state, and subsequently in a chain of 
maser transitions. 

For this inversion mechanism to operate, the pump- 
ing rate must be independent of the J quantum number. 
This is expected for collisional excitations, but also for ra- 
diative excitations if they are optically thin. The problem 
for the radiative models, of course, resides in producing 



at the same time Au=l ro- vibrational transitions opti- 
cally thin (excitation wise) and thick (de-excitation wise) . 
This is solved assuming a large velocity gradient in the 
masing region or, more generally, a particular geometry 
in which the masers arise from a thin shell. Under these 
circumstances, the radial opacity in the Ai>=l lines can 
be optically thin, while the overall average along all direc- 
tions is optically thick. In spite of this difficulty, radiative 
models offer the advantage that the vibrational excita- 
tion is much more probable than the corresponding col- 
lisional excitation. Moreover, it can naturally explain the 
thin rings of spots shown in VLBI observations, because 
the amplification is tangential, and the tight correlation 
between the SiO maser intensity and the NIR stellar ra- 
diation (Bujarrabal et al. ITM7I Pardo et al. l2T1Mjl . On 
the other hand, collisional excitation operates under quite 
general conditions, provided that the temperatures exceed 
-1000K. 

Because of the nature of the inversion mechanism, the 
physical conditions for the generation of SiO masers within 
the same vibrational state (i.e. J=l-0, J=2-l, etc.) are 
very similar regardless of the model. They occur when the 
corresponding v —* v-1 lines become optically thick, in 
a similar way for the different J levels. On the contrary, 
for masers in different vibrationally excited estates (such 
as v=l and v=2, etc.), the conditions must be naturally 
different, since it is difficult to have similar opacities for 
the Av—1 ro- vibrational transitions of levels with very 
different excitation, simply because the population of the 
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v—1 levels at ~ 1780 K will always be higher than that of 
the v—2, for the expected temperature in the inner shells 
of the envelope (3400-1700 K). This is even more stringent 
for the radiative models, since to operate, the opacity of 
the ro- vibrational transitions must be ~ 1, in order to be 
less than one in the radial direction but larger than one 
in average. 

Because of these general results, simultaneous VLBI 
observations of several maser lines, both with similar and 
different excitation energies, in the same source offer the 
opportunity to better constrain the models, by compar- 
ing their predictions with the comparison of the maps ob- 
tained for the different lines. 

As described in previous sections, we have found a dif- 
ferent result for the two sources studied. In the case of 
xCyg, the J— 2-1 line is stronger than the J=l-0, and 
the emitting regions have comparable sizes. For the v=l 
and v~2 J=l-0 lines, the v—1 maser is also stronger than 
the v=2. In IRC +10011, the differences found between the 
v=l and v=2 J=l-0 lines are smaller compared to those 
between the v=l J=l-0 and J=2-l transitions. In fact, 
the J=l-0 maser emissions are quite similar, although 
not identical, with the v=l being stronger and more ex- 
tended. On the contrary, the v=l transitions differ a lot, 
since we have obtained different angular sizes for both 
distributions, with the J=2-l one being a ~ 40% larger. 
These results for IRC+ 10011 are incompatible with the 
theoretical predictions. 

4.4. The effects of line overlaps 

One possible explanation for the different results found for 
IRC +10011 and X^yg, when comparing the 7 mm and 
3mm SiO maser emissions (see Sect. 4.1 and 4.2), may 
reside in their different chemical composition, because of 
the effects of H2O in the SiO excitation via the overlap of 
their ro- vibrational lines at 8 //m. 

It is well known that line overlap can strongly af- 
fect the pumping of maser lines because their popula- 
tion inversion is very sensitive to small changes in the 
population/de-population rates of the levels involved in 
the maser transitions. In particular, for SiO masers, it was 
already proposed by Olofsson et al. I|1981(l that the over- 
lap between the v=l, J—Q-v=2, J=l line of 28 SiO and the 
^2=0, 127_5 -z/2=l, 116,6 line of para-H20 is responsible for 
the anomalous weakness of the v=2 J=2-l SiO maser in 
O-rich stars. These authors suggested this process in view 
of the close frequencies of these two ro-vibrational tran- 
sitions, less than lkms -1 , and of the relatively strong 
v=2 J '—2-1 lines detected in S-type stars, for which lower 
H2O abundances (compared to O-rich variables) are ex- 
pected. This phenomenon was studied in more detail by 
Bujarrabal et al. I|1996l BASCC hereafter), who found 
that introducing this line overlap in the model calcula- 
tions, the observed properties of the v=2 J— 2-1 maser of 
28 SiO in both O-rich and S-type stars could be explained. 



The mechanism controlling the strength of the v—2 
3=2-1 maser via line overlap is very simple. Provided 
that the population of the v^=2 states of H2O is dom- 
inated by collisions, the opacity and excitation of the 
1/2=0, 127,5 - 1*2=1, 116,6 °f para-H 2 are larger than those 
for the v=l, J=0-v=2, 3=1 line of SiO. This results in 
an enhanced net population of the v=2, J=l SiO level, 
leading to an overthrow of the inverted population and 
the quenching off the v=2 J=2-l maser (see BASCC for 
more details). For S-type stars, however, since they are 
characterized by a [0]/[C] ratio ~ 1, the little Oxygen left 
in the gas phase after the formation of CO, should result 
in a lower H2O abundance (no H2O masers have been de- 
tected in S-type stars so far) . One must conclude that the 
properties of SiO masers in S-type stars should be easier 
to reproduce by the standard models than those of O-rich 
stars, for which the effects of overlaps could be more im- 
portant. 

This also seems to be the case in our VLBI observa- 
tions. For our M-type target, IRC +10011, for which the 
effects of the line overlap can not be neglected (as demon- 
strated by the non detection of the v=2 J=2-l maser), 
the properties of the maser spots distribution are in to- 
tal disagreement with the predictions of standard models: 
similar properties for rotational transitions in the same vi- 
brational state. For \ Cyg, in which the line overlap should 
not be as important as in M-type stars, (as demonstrated 
by the strength of its v=2 J =2-1 line), the properties of 
the maser spot distributions agree somewhat better with 
the expectations from the same standard models. 

In order to verify the influence of this line over- 
lap in the location of the masers in the envelopes of 
O-rich sources, and to test whether this effect could ex- 
plain the results obtained for IRC +10011, we have per- 
formed new calculations using the model developed by 
BASCC for SiO masers, which includes the effects of the 
1/2=0,127^-1/2=1,116,6 line of para-H 2 0. We have stud- 
ied the physical conditions under which the different SiO 
maser lines arise, for a simplified envelope model, taking 
and not taking into account the proposed line overlap. 
The characteristics of the numerical model, as well as the 
physical parameters of the envelope are the same used in 
BASCC, except for the ratio between the ro-vibrational 
and pure- rotational collisional coefficients of H2O, which 
has been assumed to be 0.02 here (a value in the range 
studied by BASCC). The results of our calculations are 
summarized in Fig. 1121 where we plot the number of maser 
photons for the v=l and 2, J=l-0 and 2-1 lines as a func- 
tion of the distance to the star and gas density, introducing 
(right panel) and not introducing (left panel) the effects 
of the overlap. 

As we can see, in the absence of overlap, the physi- 
cal conditions required for lines in the same vibrational 
state are very similar for the 7 mm and 3 mm masers for 
both v=l and v=2 transitions. In addition, the t>=2 lines 
require higher densities (and column densities) than the 
v=l masers, conditions that are very restrictive specially 
for envelopes with low mass loss rates (~ 10~ 7 Mq yr -1 )- 
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Fig. 12. Model calculations of the excitation of the SiO molecule. Number of maser photons as a function of the 
distance to the star and gas density, taking (right) and not taking into account (left) the effects of the line overlap 
between the ro-vibrational lines ^2=0,127^-^2 = 1,116.6 of para-B^O and v=l, J=0-v=2, J=l of 28 SiO (greys: v=l 
J=l-0, continuous contours: v=2 J— 1-0, dashed contours: v=l J=2— 1, dotted contours: v=2 J— 2-1). The thick grey 
lines mark the upper and lower limits of the mass loss rates in the sample of variable stars studied in BASCC, and 
the light grey dashed lines mark the mass loss rate of xCyg and IRC +10011 respectively (see Tabic [fl. 



When the overlap is included several changes occur: 

a: The v=2 J=2-l maser is quenched, except for a very 
small area in the space of physical conditions, in which 
the effect of the overlap is in fact reversed because of 
the relatively low excitation temperature of the H2O 
line (sec BASCC). 

b: The v=l J =2-1 line is barely affected by the overlap. 

c: For conditions expected in the inner shells of the enve- 
lope, especially at low mass loss rates, the v=2 J=l-0 
line requires lower densities (by a factor 2-3) to be 
masing. 

d: The conditions for the v=l J=l-0 to occur are now 
in between those for the v—1 J— 2-1 and v=2 J=l-0. 
This line is less coupled with the v—1 J— 2-1 and be- 
comes more similar to the v=2 J =1-0 line, as it is 
observed in IRC +10011. This trend is more important 
for higher values of the ratio between the ro-vibrational 
and pure-rotational collisional coefficients of H2O: for 
values of this parameter of 0.03, the v=l and v=2 
J=l-0 masing regimes become almost identical. 

Although the modeling of the envelope in this prelimi- 
nary study is very simple, the results are nevertheless very 
encouraging, since they predict more or less what we have 
observed in IRC +10011. However we believe that before 
investing more effort in this theoretical modeling more 
observational data on the relative distribution of the SiO 
maser lines, both in O-rich and S-type stars, are manda- 



tory, as well as accurate estimations of the ro-vibrational 
collisional coefficients of H2O with H2. 



5. Conclusions 

We have performed high resolution observations of masers 
of the 28 SiO molecule, in particular the rotational tran- 
sitions v=l J=l 0, v=l J=2-l, v=2 J=l-0 and v=2 
J=2 1 in the regular variables IRC +10011 and xCyg. 
The main goal of this work is to map each of these transi- 
tions so as to compare the spatial distribution of the dif- 
ferent masers. Our results allow us to study in detail the 
differences between the observational results and the pre- 
dictions of the theoretical pumping models, proving their 
validity in different types of sources. 

We have presented in this paper the first VLBI maps 
for an S-type star, xCyg: in particular, the 28 SiO rota- 
tional lines J =1-0 and J=2-l in the v=l and v=2 vibra- 
tional states. For the O-rich star, IRC +10011, we have 
mapped the 28 SiO v=l J=l-0 and J=2-l, and the v=2 
J =1—0 maser emissions. 

In IRC +10011, all the 28 SiO maser transitions were 
detected and mapped except the V=2 .7=2 1 line. When 
comparing the two 7 mm maser emissions, we found that 
the v=2 is located in an inner region of the envelope, its 
angular extent being 10% smaller than that of the v=2. 
On the contrary, for the v=l 7=1-0 and 7=2-1 maser 
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lines, the latter emission arises in a layer ~50% farther 
away from the star. 

For xCyg, the v=l J=l-0 and J=2-l have compara- 
ble sizes, although their spatial distributions are different. 
On the other hand, the v=2 lines, J=l-0 and J =2-1, seem 
to originate in a similar region probably more compact and 
with more restrictive characteristics in the envelope. 

Comparing our observational results with the different 
theoretical pumping mechanisms we obtain the following: 

- Standard models, collisional or radiative, have se- 
rious problems in explaining what is observed in 
IRC+lOOff: the v=l and v=2 distributions being 
more similar than those of the v=l J=l-0 and J=2-l 
transitions. This observational result contradicts what 
is predicted by all the theoretical studies, where the 
distributions of the J=l-0 and J =2-1 emissions ap- 
pear to be similar. 

- However, in \ Cyg, although the distributions of the 
v=l J=l-0 and J =2—1 transitions are different, as in 
IRC+lOOff, the sizes of the masing regions are com- 
parable. For the v=2 J=l-0 and J =2-1 lines, it is 
difficult to establish the differences between the two 
emitting regions since, in the first transition only one 
feature was found. 

- With respect to the J =2-1 and J =1-0 maser lines, 
a special effort must be made to study whether the 
observed relative spatial distribution of these emissions 
is similar to IRC +10011 or \ Cyg- 

- The fact that almost none of the spots are coincident 
in all the transitions reveals that there may be addi- 
tional maser competition effects that are not theoret- 
ically understood. We note that, since a considerable 
percentage of the emission is lost, the lack of coinci- 
dence could refer only to the brightest features. 

To explain the discrepancies between theory and ob- 
servations, we have proposed that the line overlap between 
ro- vibrational H2O and SiO lines is a process that strongly 
affects the conditions under which the masers occur. This 
effect is stronger in O-rich variables, such as IRC +10011, 
since they present a higher abundance of H2O in the enve- 
lope than S-type variables, such as \ Cyg. The results are 
very promising: when the overlap is important, the v=l 
and v=2 J=l- lines are strongly coupled and a signif- 
icant difference appears between the J=l-0 and J=2 1 
maser lines, as is observed in IRC +10011. 

We therefore conclude that a more extended study, 
from observational and theoretical points of view, will help 
to understand if the results for the SiO maser emissions 
obtained in IRC +10011 and x^yg are produced system- 
atically in other LPVs and to constrain both the pumping 
mechanisms and the importance of the line overlaps. 
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